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Fractional Cauchy problems replace the usual first-order time 
derivative by a fractional derivative. This paper develops classical 
solutions and stochastic analogues for fractional Cauchy problems 
s ! ■ in a bounded domain D C R d with Dirichlet boundary conditions. 

Stochastic solutions are constructed via an inverse stable subordi- 
nator whose scaling index corresponds to the order of the fractional 
time derivative. Dirichlet problems corresponding to iterated Brow- 
nian motion in a bounded domain are then solved by establishing a 
correspondence with the case of a half-derivative in time. 

1. Introduction. In this paper, we extend the approach of Meerschaert 
and Schefner [23] and Meerschaert et al. [24] to fractional Cauchy problems 
£S) . on bounded domains. Our methods involve eigenfunction expansions, killed 

Markov processes and inverse stable subordinators. In a recent related paper 
[7] , we establish a connection between fractional Cauchy problems with index 
\q ', j3 = 1/2 on an unbounded domain, and iterated Brownian motion (IBM), 

' defined as Zt = B(\Yt\), where B is a Brownian motion with values in M. d 

£\\ ■ and Y is an independent one-dimensional Brownian motion. Since IBM is 

also the stochastic solution to a Cauchy problem involving a fourth-order 
derivative in space [2, 14], that paper also establishes a connection between 
certain higher-order Cauchy problems and their time-fractional analogues. 
More generally, Baeumer, Meerschaert and Nane [7] shows a connection 
^ ■ between fractional Cauchy problems with (3 = 1/2 and higher-order Cauchy 

problems that involve the square of the generator. In the present paper, we 
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clarify that relationship by identifying the boundary conditions that make 
these two formulations have the same unique solutions on bounded domains; 
see Corollary 4.5. 

A celebrated paper of Einstein [15] established a mathematical link be- 
tween random walks, the diffusion equation and Brownian motion. The scal- 
ing limits of a simple random walk with mean zero, finite variance jumps 
yields a Brownian motion. The probability densities of the Brownian mo- 
tion variables solve a diffusion equation, and hence we refer to the Brownian 
motion as the stochastic solution to the diffusion equation. The interplay 
between stochastic processes and partial differential equations has many 
practical and theoretical applications. For example, it has lead to efficient 
particle tracking simulation methods for solving complicated boundary value 
problems [18]. The diffusion equation is the most familiar Cauchy problem. 
The general abstract Cauchy problem is dtu = Lu, where u(t) takes values 
in a Banach space and L is the generator of a continuous semigroup on that 
space [4]. If L generates a Markov process, then we call this Markov process 
a stochastic solution to the Cauchy problem dtu = Lu, since its probability 
densities (or distributions) solve the Cauchy problem. This point of view 
has proven useful, for instance, in the modern theory of fractional calculus, 
since fractional derivatives are generators of certain (a-stable) stochastic 
processes [23]. 

Fractional derivatives are almost as old as their more familiar integer- 
order counterparts [26, 30]. Fractional diffusion equations have recently been 
applied to problems in physics, finance, hydrology and many other areas 
[17, 20, 25, 32]. Fractional space derivatives are used to model anomalous 
diffusion or dispersion, where a particle plume spreads at a rate inconsistent 
with the classical Brownian motion model, and the plume may be asym- 
metric. When a fractional derivative replaces the second derivative in a 
diffusion or dispersion model, it leads to enhanced diffusion (also called 
superdiffusion). Fractional time derivatives are connected with anomalous 
subdiffusion, where a cloud of particles spreads more slowly than a classical 
diffusion. Fractional Cauchy problems replace the integer time derivative by 
its fractional counterpart: d^u = Lu. Here, d^g(t) indicates the Caputo frac- 
tional derivative in time, the inverse Laplace transform of s^g(s) — s^~ 1 g(0), 
where g(s) = e~ st g(t) dt is the usual Laplace transform [11]. Nigmatullin 
[27] gave a physical derivation of the fractional Cauchy problem, when L is 
the generator of some continuous Markov process {Y(t)} started at x = 0. 
The mathematical study of fractional Cauchy problems was initiated by 
[19, 20, 33]. The existence and uniqueness of solutions was proved in [19, 20]. 
Fractional Cauchy problems were also invented independently by Zaslavsky 
[35] as a model for Hamiltonian chaos. 

Stochastic solutions of fractional Cauchy problems are subordinated pro- 
cesses. If X(t) is a stochastic solution to the Cauchy problem dtu = Au, 
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then under certain technical conditions, the subordinate process X{Et) is 
a stochastic solution to the fractional Cauchy problem d^u = Au; see [5]. 
Here, Et is the inverse or hitting time process to a stable subordinator Dt 
with index (3 S (0,1). That is, Et = inf{a; > 0:D X > t}, and Dt is a Levy 
process (continuous in probability with independent, stationary increments) 
whose smooth probability density gp(t) has Laplace transform e~ s = gp(s); 
see [31]. Just as Brownian motion is a scaling limit of a simple random 
walk, the stochastic solution to certain fractional Cauchy problems are scal- 
ing limits of continuous time random walks, in which the i.i.d. jumps are 
separated by i.i.d. waiting times [24]. Fractional time derivatives arise from 
power law waiting times, where the probability of waiting longer than time 
t > falls off like t~@ for t large [23]. This is related to the fact that fractional 
derivatives are nonlocal operators defined by convolution with a power law 
[5]. 

2. Preliminaries. Let D be a bounded domain in W 1 . We define the 
following spaces of functions: 

C(D) = {u : D — > K : u is continuous}, 

C(D) = {u:D ^ IR : u is uniformly continuous}, 
C k {D) = {u:D^R :u is fc-times continuously differentiable}, 
C k (D) = {u£C k (D):D~< u is uniformly continuous for all |7| < k}. 

Thus, if u € C k (D), then D^u continuously extends to D for each multi- 
index 7 with |7| < k. 

We define the spaces of functions C°°(D) = f]^ =1 C k {D) and C°°(£>) = 
f]T=iC k (D). 

Also, let C k ' a (D) [C k ' a (D)] be the subspace of C k {D) [C k (D)] that con- 
sists of functions whose kth order partial derivatives are uniformly Holder 
continuous with exponent a in D. For simplicity, we will write 

C°^(D) = C a (D), C Q ' a {D) = C a (D) 

with the understanding that < a < 1 whenever this notation is used, unless 
otherwise stated. 

We use C c {D),C k {D),C k > a {D) to denote those functions in C(D), C k {D), 
C k,a (D) with compact support. 

A subset D of M. d is an /-dimensional manifold with boundary if every 
point of D possesses a neighborhood diffeomorphic to an open set in the 
space H l , which is the upper half space in M 1 . Such a diffeomorphism is 
called a local parametrization of D. The boundary of D, denoted by dD, 
consists of those points that belong to the image of the boundary of H l 
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under some local parametrization. If the diffeomorphism and its inverse are 
(jk,a f unc tions, then we write dD G C k,a . 

We say that D satisfies an exterior cone condition at a fixed point xq G dD 
if there exists a finite right circular cone V = V XQ with vertex xq such that 
D H V XQ = xo , and a uniform exterior cone condition if D satisfies an exterior 
cone condition at every point xq G dD and the cones V xo are all congruent 
to some fixed cone V. 

Let {X t } be a Brownian motion in M rf and A C M d . Let Ta = inf{i > 
0: G A} be the first hitting time of the set A. We say that a point y is 
regular for a set A if iy[Tk = 0] = 1. Note that a point y is regular for a set 
A provided, starting at y, the process does not go a positive length of time 
before hitting A. 

The right condition for the existence of the solution to the Dirichlet prob- 
lem turns out to be that every point of dD is regular for D c (cf. [8], Section 
HI)- 

If a domain satisfies a uniform exterior cone condition, then every point 
of dD is regular for D c . 

The Sobolev space W k,p (D) consists of all locally integrable functions 
u:D — ► M. such that for each multi-index 7 with I7I < k, D^u exists in the 
weak sense and belongs to L P (D) with the following norm: 



W k 'P(D) — I \\J\\p,D 



i 7 i<fe / 



i/p 



We denote by Wq' p (D) the closure of C k {D) in W k 'P(D). If p = 2, we 
usually write 



H k {D) = W k ' 2 (D), H$(JD) = WS'*(D) (k = 0, 1, . . .)• 

Let dD G C . Then at each point x G dD there exists a unique outward 
pointing unit vector 

6(x) = (e 1 (x),...,0 d (x)). 

Let u G C 1 (-D), the set of functions which have continuous extension of 
the first derivative up to the boundary. Let 

„ du _ _ 
D 6 u = — = e-Vu, 

denote the directional derivative, where Vm is the gradient vector of u. 

Now we recall Green's first and second identities (see, e.g., [16], Section 
2.4). Let u,veC 2 (D)r\C 1 {D). Then 
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/ — — vdx= [ uvOids— [ u— — dx 
Id oxi Jqd Jd oxi 



(integration by parts formula), 

dv 



[ Vv-Vudx = — [ uAvdx+ [ —uds (Green's first identity), 
Jd Jd JdD oO 

/ [uAv — vAu] dx = 
Jd Ja 



ID JD JdD 

ds (Green's second identity). 



dv du 

u m~ v m 



ID JdD 

A uniformly elliptic operator of divergence form is defined on C 2 functions 

by 

.4-"^ dxj 

1,3=1 ■> 
with aij(x) = aji(x) and, for some A > 0, 

n n n 

(2.2) A^ % 2 < ^ aij (x)y iyj < A" 1 £ yf Vy G R d . 

i=l i,j=l i=l 

The operator L acts on the Hilbert space L 2 (D). We define the initial 
domain Cq°(D) of the operator as follows. We say that / is in Cq°(D), if 
/ G C°°{D) and f(x) = for all x G dD. This condition incorporates the 
notion of Dirichlet boundary conditions. 

From [13], Corollary 6.1, we have that the associated quadratic form 



is closable on the domain Cq°(D) and the domain of the closure is inde- 
pendent of the particular coefficients (««) chosen. In particular, for /, g G 
C^°(D) by integration by parts 



g(x)Lf(x)dx = Q(f,g)= / f(x)Lg(x)dx, 
d Jd 

which shows that L is symmetric. 

From now on, we will use the symbol Lo if we particularly want to empha- 
size the choice of Dirichlet boundary conditions, to refer to the self-adjoint 
operator associated with the closure of the quadratic form above by the use 
of [13], Theorem 4.4.5. Thus, Ld is the Friedrichs extension of the operator 
defined initially on Cq°(D). 
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If the coefficients dij(x) are smooth \aij{x) 6 C 1 (D)], then Lj^u takes the 
form 

Ldu= £ Mx) a^ + ^ffaatW) *i 



dxj dx~ ^—f \ f-f &c,- / <9x 

t,J=l ' J «=1 Yj=l J 



d 



d 2 u , . . du 



If Xt is a solution to 



dxndxj ' ^-rf tv dxi 

i,3=l J t=l 



dX t = a(X t )dW t + b{X t )dt, X = x , 

where a is a d x (i matrix, and Wt is a Brownian motion, then Xt is associated 
with the operator Ld with a = aa T (see Chapters 1 and 5 of Bass [8]). Define 
the first exit time as td(X) = mf{t >0:X t <£ D}. The semigroup defined by 
T(t)f(x) = E x [f '{Xt)l (td(X)) > £)] has generator Ld, which follows by an 
application of the Ito formula. 

Let D be a bounded domain in M. d . Suppose L is a uniformly elliptic 
operator of divergence form with Dirichlet boundary conditions on D, and 
that there exists a constant A such that for all x £ D, 

d 

(2.3) |ay(aO|<A. 

Let Tb(i) be the corresponding semigroup. Then Tb(i) is an ultracontractive 
semigroup (even an intrinsically ultracontractive semigroup); see Corollary 
3.2.8, Theorems 2.1.4 and 4.2.4 and Note 4.6.10 in [12]. Every ultracontrac- 
tive semigroup has a kernel for the killed semigroup on a bounded domain 
which can be represented as a series expansion of the eigenvalues and the 
eigenfunctions of Ld (cf. [12], Theorems 2.1.4 and 2.3.6 and [16], Theorems 
8.37 and 8.38): There exist eigenvalues < fii < fi2 < ^3 < ■ • ■ , such that 
li n — > oo, as n — > oo, with the corresponding complete orthonormal set (in 
Hq ) of eigenfunctions tp n of the operator L d satisfying 



(2.4) L D ip n (x) = -fJLn1pn(x), X G D : 1pn\dD = 0. 

In this case, 

oo 

p D (t,x,y) = e -,tet i(i)^(i/) 

n=l 

is the heat kernel of the killed semigroup Tp . The series converges absolutely 
and uniformly on [to, oo) x D x D for all to > 0. 
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Denote the Laplace transform by 

POO 

u(s,x) = / e~ st u(t,x) dt. 



o 



Since we are working on a bounded domain, the Fourier transform methods 
in [24] are not useful. Instead we will employ Hilbert space methods. Hence, 
given a complete orthonormal basis {ij} n (x)} on Hq, we will call 



u(t,n)= / ip n {x)u(t,x) dx, 
Jd 

u(s,n) = / ip n (x) / e~ st u(t,x) dtdx 



D 

Ms) I e" st - 

D 

t/j n (x)u(s, x) dx 



D 

the ip n and V'n-Laplace transforms, respectively. Since {tp n } is a complete 
orthonormal basis for Hq, we can invert the ^-transform 

u(t,x) =^2u(t,n)ip n (x) 

n 

for any t > 0, where the sum converges in the 1? sense (e.g., see [29], Propo- 
sition 10.8.27). 

Suppose D satisfies a uniform exterior cone condition. Let {Xt} be a 
Markov process in R d with generator Lp, and / be continuous on D. Then 
the semigroup 

T D (t)f(x) = E x [f(X t )I(t < t d (X))} 
(2-5) = / p D (t,x,y)f(y)dy 



D 

oo 



= ^e-"»V„(x)/(n) 

n=l 

solves the Dirichlet initial-boundary value problem in D: 

du{t,x) =LdU ^ x ^ x eD,t>0, 

u(t,x) =0, x E dD, 
u(0,x)=f(x), x G D. 

Remark 2.1. The eigenfunctions belong to L°°(D) n C a (D) for some 
a > 0, by [16], Theorems 8.15 and 8.24. If D satisfies a uniform exterior cone 
condition all the eigenfunctions belong to C a (D) by [16], Theorem 8.29. If 
Oij G C a {D) and <9D G C 1,ct , then all the eigenfunctions belong to C 1,a {D) 
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by [16], Corollary 8.36. If ctjj G C°°(D) then each eigenfunction of L is in 
C°°(D) by [16], Corollary 8.11. If ay G C°°(Z)) and dD eC°°, then each 
eigenfunction of L is in C°°(D) by [16], Theorem 8.13. 

([16], Corollary 8.36). Let D be bounded and every point of dD be regular 
for D c . In the case = A^, the corresponding Markov process is a killed 
Brownian motion. We denote the eigenvalues and the eigenfunctions of Ad 
by {A n , 4> n } 1 , where (f> n G C°°(D). The corresponding heat kernel is given 
by 

oo 

PD(t,x,y) = J2e~ Xnt M%)My)- 

n=l 

The series converges absolutely and uniformly on [to, oo) x D x D for all 
to > 0. In this case, the semigroup given by 

T D (t)f(x) = E x [f(X t )I(t < t d (X))] 

r oo 

= / p D (t,x,y)f(y)dy = Y,e~ Xnt Mx)f(n) 

solves the heat equation in D with Dirichlet boundary conditions: 

du(t,x) . . . „ 

— ^ = Au(t,s), x€D,t>0, 

u(t,x) =0, x G dD, 

u(0,x) = /(x), x G D. 

The Caputo fractional derivative of order < /3 < 1 , defined by 

, . 9^n(t,x) 1 [ l du(s,x) ds 

(2.6) 



9^ r(l-/3)7o 9s (t-s)/ 3 ' 

was invented to properly handle initial values [11]. Its Laplace transform 
sPu(s,x) — s^ -1 u(0,x) incorporates the initial value in the same way as the 
first derivative. The Caputo derivative has been widely used to solve ordinary 
differential equations that involve a fractional time derivative [17, 28]. 

Inverse stable subordinators arise in [23, 24] as scaling limits of contin- 
uous time random walks. Let S(n) = Y\ + • • • + Y n a sum of i.i.d. random 
variables with EY n = and EY% < oo. The scaling limit c -1 / 2 : S([ct]) =^ B(t) 
as c — > oo is a Brownian motion, normal with mean zero and variance pro- 
portional to t. Consider Y n to be the random jumps of a particle. If we 
impose a random waiting time T n before the nth jump Y n , then the position 
of the particle at time T n = J\ + • • • + J n is given by S(n). The number of 
jumps by time t > is Nt = max{n:T„ < t}, so the position of the parti- 
cle at time t > is S(N t ), a subordinated process. If P(J n >t) = t~ L(t) 
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for some < (3 < 1, where Lit) is slowly varying, then the scaling limit 
c^'^TifM => Dt is a strictly increasing stable Levy motion with index j3, 
sometimes called a stable subordinator. The jump times T n and the num- 
ber of jumps Nt are inverses {Nt > x} = {T(\x}) < t}, where \x] is the 
smallest integer greater than or equal to x. It follows that the scaling 
limits are also inverses c~@N c t Et, where Et = inf{a;:-D x > t}, so that 
{Et < x} = {D x > t}. Since N c t ~ c^Et, the scaling limit of the particle 
locations is c _/3/2 5(iV[ ci ]) « (c^)" 1 / 2 S(cP E t ) « B(E t ), a Brownian motion 
subordinated to the inverse or hitting time (or first passage time) process of 
the stable subordinator Dt- The random variable Dt has a smooth density. 
For properly scaled waiting times, the density of Dt has Laplace transform 
g-ts* 9 £ Qr an y ^ y an j jj^ j g identically distributed with t l ^D\. Writing 
gp(u) for the density of D\, it follows that D t has density t~ 1 /^ gp{t~ l /^u) 
for any t > 0. Using the inverse relation -P(£t < x) = P(D X > t) and taking 
derivatives, it follows that Et has density 

(2.7) f t (x)=tp- 1 x- 1 - 1/P 9p(tx- 1 /P), 

whose 1 1 — ^ s Laplace transform s@~ 1 e~ xs ^ can also be derived from the equa- 
tion 

ft(x) = ^P(D X > t) = ^j™x- 1 !( 3 g f} {x- 1 'f 3 u)du 
by taking Laplace transforms on both sides. 

3. Fractional diffusion in bounded domains. Fractional Cauchy problems 
replace the usual first-order time derivative with its fractional analogue. In 
this section, we prove classical (strong) solutions to fractional Cauchy prob- 
lems on bounded domains in M. d . We also give an explicit solution formula, 
based on the solution of the corresponding Cauchy problem. Our methods 
are inspired by the approach in [24], where Laplace transforms are used 
to handle the fractional time derivative, and spatial derivative operators (or 
more generally, pseudo-differential operators) are treated with Fourier trans- 
forms. In the present paper, we use an eigenfunction expansion in place of 
Fourier transforms, since we are operating on a bounded domain. Our first 
result, Theorem 3.1, is focused on a fractional diffusion d^u = Au, and we 
lay out all the details of the argument in the most familiar setting. Then, in 
Theorem 3.6, we use a separation of variables to handle arbitrary generators: 
dfU = Lu. Along the way, we explicate the stochastic solutions in terms of 
killed Markov processes. 

Let (3 G (0, 1), L>oo = (0, oo)x D and define 
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n A (D c 



Of 



u(t, x) 



d d & 

<g(x)tP-\ g£L°°(D), t > o|. 



We will write u G C k (D) to mean that for each fixed t > 0, u(t, •) G C k (D), 
and w G (AxO to mean that it G C k (D OQ ) and is bounded. 

Theorem 3.1. Let < a < 1. Ze£ D be a bounded domain with dD G 
C 1,Q ; and Tu(t) be the killed semigroup of Brownian motion {Xt} in D. 
Let Ef be the process inverse to a stable subordinator of index (3 G (0, 1) 
independent of {X t }. Let f G D{A D ) n C l (D) n C 2 (D) for which the eigen- 
function expansion ( of Af ) with respect to the complete orthonormal basis 
{4> n :nGN} converges uniformly and absolutely. Then the unique (classical) 
solution of 



(3.1) 



u g WaPoo) n c 6 (Ax>) n c\D), 

-gjpu(t,x) = Au(t,x), x£D,t>0, 



u(t,x) = 0, xedD,t>0, 
u(0, x) = f(x), x G D 

is given by 

u(t,x) = E x [f(X{E t ))L{T D {X) > Et)] 

= ~0Jo T ^f( x )9^~ 1/P T 1/P ~ 1 ^ = J T D ((t/lf)f(x)g p (l)dl. 

Proof. Assume that u(t,x) solves (3.1). Using Green's second identity, 
we obtain 



/ [uAcj) n - 4> n Au] dx = 
Jd Jd 



dD 
*1< 



!/■■ 



de 



du 
86 



ds = 0, 



since u\qd = = 4> n \dD and u,cf) n G C (D). Hence, the ^-transform of Au 
is 



4> n (x)Au(t,x) dx = / u(t, x)A(j) n (x) dx 
d Jd 



-A n / u(t,x)(j) n (x)dx = -\ n u(t,n), 
Jd 
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as (f) n is the eigenfunction of the Laplacian corresponding to eigenvalue A n . 

Recall the definition (2.6) of the fractional derivative, and use the fact 
that | du ^ x > | is bounded by g(x)t^~ 1 where g G L°°(D). It follows that the 
fractional derivative commutes with (^-transform, and hence, we can take 
4> n -transforms in (3.1) to obtain 

(3.2) -^u(t,n) = -\ n u(t,n). 

Since u is uniformly continuous on C([0, e]) x D) and hence uniformly bounded 
on [0, e] x D, we have by dominated convergence that lim t ^o J D u(t,x) 
4> n (x) dx = f(n). Hence, u(0,n) = f(n). A similar argument shows that t \— > 
u(t, n) is a continuous function of t G [0, oo) for every n. Then, taking Laplace 
transforms on both sides of (3.2), we get 

(3.3) sPu{8, n) - s /3 ~ 1 n(0, n) = -X n u{s, n) 
which leads to 

(3.4) «.,») 

By inverting the above Laplace transform, we obtain 

u(t,n)=f(n)E f3 {-\ n tf 3 ) 
in terms of the Mittag-Lefner function defined by 

oo z k 

see, for example, [22]. Inverting the ^-transform, we get an L 2 -convergent 
solution of equation (3.1) as 



(3.5) u(t,x) = f{n)Mx)Ep{-\ n t^ 

n=l 



for each t > 0. In order to complete the proof, it will suffice to show that the 
series (3.5) converges pointwise, and satisfies all the conditions in (3.1). 

Step 1. We begin showing that (3.5) convergence uniformly in t G [0, oo) 
in the L 2 sense. Define the sequence of functions 

N 

(3.6) u N (t, x) = J2 /(n)<M*)£/3(-A„A 

71=1 

We have, by [21], equation (13), 

(3.7) 0<^(-A„^)<c/(l + A n t /3 ). 
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Since / 6 L 2 (D), we can write f(x) = Y^=i f{ n )<fin{x), and then the Par- 
seval identity yields 



£(/») 2 = II/III,d<°°- 

71=1 

Then, given e > 0, we can choose no(e) such that 

oo 

(3.8) £ (/») 2 < e - 

n=no(e) 

For N > M > n (e) and i > 0, 

AT 



\u N (t,x) - u M (t,x) 



\2,D 



(3.9) 



< 



< 



n=M 



2,D 



(1 + A no ^) 



£ (/») 2 

n=no (e) 



< c 2 E (/») 2 <^- 

n=re (e) 

Thus, the series (3.5) converges in L 2 (D), uniformly in t > 0. 

Step 2. Next we show that the initial function is defined as the L 2 limit 
of u(t, x) as t — > 0, that is, we show that t — > u(-,t) 6 C((0, oo); L 2 (D)) and 
u(-,t) takes the initial datum / in the sense of L 2 (D), that is, 

N-,t)-/||2,r>->0 ast^O. 

Note Ep(—z) is completely monotone for z > (cf. Agrawal [1]). In partic- 
ular, Ep(—\ n tP) is a nonincreasing function of A n . Hence, we have 



t(t,x) - /(x) = E /W(^(-A n ^) " l)0n(s). 

Fix e E (0, 1) and choose no = no(e) as in (3.8). Then, 

oo 

■,*) - f\\l,D = £(/») 2 CM-^) - 1) 2 



71=1 

no (e) 

< £ (/») 2 (%(-A^) - if 



71=1 



+ e (/») 2 (%(-A7i^)-ir 

n=7io(e)+l 

<(l-^(-A no ^)) 2 ||/|| 2 ,z) + e 
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and now the claim follows, since Ep{—\ n t^) — ► 1, as t — > 0. The continuity 
of t \-+u(t, •) in L 2 (D) at every point t £ (0, oo) can be proved in a similar 
fashion. 

Step 3. We now obtain a decay estimate for the solution u(t,x). Using 
Parseval's identity, the fact that A n is increasing in n, and the fact that 
Ep(—\ n tP) is nonincr easing for n > 1, we get 

|K^)||2,D<%(-Al^)||/|| 2 ,D. 

Step 4. We show uniform and absolute convergence of the series (3.5) 
defining u(t,x), which will imply that the classical solution to (3.1) is given 
by this equation. We do this by showing that (3.6) is a Cauchy sequence in 
L°°(D) uniformly in t > 0. 

Applying the Green's second identity we see that 



A/(n) 



D 



Af(x)(j) n (x) dx = (-A n )/(n). 



Hence the series defining Af is J2^=i — A„0„(x)/(n). So, this series is abso- 
lutely and uniformly convergent. Let e > 0. Since \ n — ► oo as n — > oo, there 
exists an no(e) so that for all x € D, 

oo oo 

(3.10) J2 l/»ll<M*)|< E A n |/"(n)||0 n (x)|<6. 

n=no(e) n=no(e) 

This is possible since we have assumed that Af has an expansion which is 
uniformly and absolutely convergent in L°°{D). We will freely use the fact 
that the series defining / also converges absolutely and uniformly. 
For N>M> n (e) and t > and x G D, 

N 

E </'n(x)/(n)^(-A n t /3 ) 



|itjv(t,x) - u M (t,x)\ 



(3.11) 



n=A*f 



< 



n T x iflN E l^n(^)||/(«)| <ce. 

,1 + An ^) n=M 



This shows that the sequence UN(t,x) is a Cauchy sequence in L°°(D) and 
so has a limit in L°°(D). Hence, the series in (3.5) is absolutely and uni- 
formly convergent. It follows that u(t,x) satisfies the boundary conditions 
in equation (3.1). 

Step 5. Next we show that the fractional time derivative 3f and the 
Laplacian A can be applied term by term in (3.5). Note that 



E }{n)Ep{-\ n t p )A^ n {x) 



n=l 



Y^f(n)Ep{-\ n tP)\ n (i> n {x) 



n=l 
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<$> n (z)||/»| r 



n=l 



<t-^l^)||/(n)|<oo, 



n=l 



where the last inequality follows from the fact that the eigenfunction expan- 
sion of / converges absolutely and uniformly. Then the series 

oo 

/(n)%(-A n ^)A^ n (x) 

n=l 

is absolutely convergent in L°°(D) uniformly in (io,oo) for any to > 0. The 
function Ep{—Xt^) is an eigenfunction of the Caputo fractional derivative, 

with d^E^-Xt 13 ) = -XEpi-Xt 13 ) for any A > 0, which follows by taking the 
Laplace transform of both sides of this equation. Hence 

oo 0/3 oo 

£ f(n)Mx) W E p (-X n t^ = ]T J{n)E p (-\ n tP)^ n (x) 

71=1 71=1 

since the two series are equal term by term, and since the series on the right 
converges. 

Now it is easy to check that the fractional time derivative and Laplacian 
can be applied term by term in (3.5) to give 

-,u{t, x) — Au(t, x) 



an 3 



n=l 



n 



'OH 3 ' 



(f )n {x)^E /3 (-X n t (3 ) - Ep{-X n t p )A<j) n (x) 



so that the PDE in (3.1) is satisfied. Thus, we conclude that u defined by 
(3.5) is a classical (strong) solution to (3.1). 
Using [21], equation (17), 

dE p (-X n tP) 



dt 



we get 



du(x,t) 



dt 



oo 

<£i/(«)i 

n=l 



<c Kt [ \ <cX n tP-\ 

~ i + x n tP- 

dEp{-X n tP) 



dt 



\4>n(x)\ 



< dP- 1 £ X n \f(n)\\Mx)\ = t^gix). 



n=l 



Since A/ has absolutely and uniformly convergent series expansion with 
respect to {4> n : n G N}, we have g(x) £ L°°(D). 
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Now from the results obtained above it follows that u G r HA{D o)^Cb{D 00 ). 

We next show that u S C 1 (L>); this follows from the bounds in [16], The- 
orem 8.33, and the absolute and uniform convergence of the series defining 
/• 

(3.13) \4>n\l,a;D < C(l + A n ) SUp \(p n {x)\ , 

D 

where C = C(d,X,A,dD) and A is the constant in the definition of uniform 
ellipticity of L in (2.2) and A is the bound in (2.3). Here 

k 

\u\k,a;D = SUp [-D 7 ti] Q ,D + ^ sup sup|D 7 u|, k = 0, 1, 2, . . . , 

M=fc j=o\j\=j D 

and 

[D^uI^d = sup 1 _ 

are norms on C k ' a (D). Hence 

oo 

H;t)\ ha . D <CV |/»|^(-A„^)(1 + A n ) sup \M*)\ 
n=l D 

00 1 I \ 

<c£|/>)|-^-^su P |0 n (*)| 

oo oo 

< Ct-P £ sup 1 n (x) 1 1 /(n) | + C ]T sup | <An (x) 1 1 / (n) | < oo . 

n=l D n=l D 

Step 6. Next we obtain the stochastic solution to (3.1), by inverting the 
^ n -Laplace transform, using a method similar to [24]. Using the fact that 
for any A > 0, 

OO 1 

e~ xl dl = \, 
A 



we get 



u(s, n) - 



(3.14) = J{n)s p - 1 r e~^ +K)l dl 

Jo 

' — ft 

e~ Xnl f(n) S P- l e- ls ^ dl. 
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The </> n -transform of the killed semigroup Tr>{l)f(x) = J2m=i e ~ Xrnl 4>m{ x ) x 
/ (m) from (2.5) is found as follows. Since {4> n , n £ N} is a complete orthonor- 
mal basis of L 2 (D), we get 



-l\ n 



D 



[T D {l)f]{n)= / Mx)T D (l)f(x)dx 
Jd 

By [24], equation (8), we have 

3 te-tgpitrwy-v+w dt, 



f{n){^ n {x)fdx = e~ lx -f(n 



1 



where gp{t) is the smooth density of the stable subordinator with Laplace 
transform gp{s) = / °° e~ st gp(t) dt = e~ sli . Combining the last two results, 
we get 

"l r c 



u(s, n) 



(3.15) 



[T D {l)f]{n) 



te- 8t gp{tt-W)rQ- +1 /® dt 



-St 



di 

dt. 



Invert the above Laplace transform to obtain 



u(t, n) 



W^W}(n)^(ti-^)r^^ di. 



Inverting the ^-transform is equal to multiplying both sides of the above 
equation by <j) n and then summing up from n = 1 to oo. 

We have uniqueness of this inverse in L 2 (D) and u(t, •) £ L 2 (D) for each 
fixed t>0. Inverting the ^-transform, we obtain 



t 



u(t,x) = £ Mx) J o [T D (/)/](n)-^(t/- 1 //3)/-(i+i//3) dL 



Using 
(3.16) 
we get 



n=l 



- x ^gp{l)dl = Ep{-\ n t^), 



°° POO + 

c(t,x) = £ Mx) / [T D (/)/](n)- 5/3 (t/- 1 ^)/-( 1+1 /' 3 ) dl 

n=l J ° P 

oo 



n=l 



After showing the absolute and uniform convergence of the series defin- 
ing u, we can use a Fubini-Tonelli type argument to interchange order of 
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summation and integration in the following to get stochastic representation 
of the solution as 

00 POO f 

£1 J ° P 



(3.17) 



Ln=l 



-l\ n 



Igpia-wy-wn di 







1 + 

Tn^m-gpitrV^l-^/Vdl 



= E x [f(X(E t ))I{T D {X)>E t )]. 

The last equality follows from a simple conditioning argument, and the den- 
sity formula (2.7) for the inverse stable subordinator. 

Step 7. Finally we prove uniqueness. Let u\, ui be two solutions of (3.1) 
with initial data u(0,x) = f(x) and Dirichlet boundary condition u(t,x) = 
for x G dD. Then U = u\ — ui is a solution of (3.1) with zero initial data 
and zero boundary value. Taking n -transform on both sides of (3.1) we get 

- 

—U( t ,n) = -\ n U(t,n), C/(0,n)=0, 

and then U(t,n) = for all t > and all n > 1. This implies that U(t,x) = 
in the sense of L 2 functions, since {<p n : n > 1} forms a complete orthonormal 
basis for L 2 (D). Hence, U (t, x) = for all t > and almost all x G D. Since U 
is a continuous function on D, we have U(t, x) = for all (t, x) E [0, 00) x D, 
thereby proving uniqueness. □ 

Corollary 3.2. The solution in Theorem (3.1) has also the following 
representation: 

u(t,x) = E x [f(X(E t ))I( TD (X) > E t )] = E x [f(X(E t ))I{T D {X(E)) > t)}. 

Proof. Given a continuous stochastic process Xt on R , and an interval 
/ C [0, 00), we denote X(I) = {X u :u 6 /}. Since the domain D is open and 
X u is continuous, it follows that 

{r D (X)>t} = {X([0,t])cD}. 

To see this, first note that X u G D for all < u < t implies that Xt G D. 
Choose an open set B C D such that Xt£B. Since X is continuous, the set 
X~ l (B) = {u : X u G B} is open, and it contains t. Then it also contains t + e 
for any e > sufficiently small, so that X u G D for all < u < t + e, and hence 
t d (X) > t. This shows that {X([0,t]) C D} C {t d (X) > t}. Conversely, if 



18 



M. M. MEERSCHAERT, E. NANE AND P. VELLAISAMY 



t d (X) > t then X u G D for all < u < t, so that {t d (X) > t} C {X([0,i]) C 
Z)}, and hence these two sets are equal. 

Next note that, since Et is continuous and monotone nondecreasing, we 
have E([Q,t]) = [0,E t ]. To see this, first note that, since E t is nondecreas- 
ing with E = 0, u G [0,t] implies E u G [0,E t ], so that E([0,t]) C [0,E t ]. 
Conversely, if y G [0,E t ] then suppose y ^ -EQO,^). Then (£^[0,i]) is not con- 
nected, which contradicts the fact that the continuous image of a connected 
set is connected. Hence [0,i?t] C E'QOji]), so these two sets are equal. 

Finally, we observe that 



which completes the proof. □ 

Remark 3.3. If we subordinate Brownian motion X by a stable subor- 
dinator Dt in the sense of Bochner, then the conclusions of Corollary 3.2 do 
not hold since a stable subordinator does not have continuous sample paths; 
see, for example, Song and Vondracek [34]. Subordination by the inverse 
stable subordinator has the property that killing the process X and then 
subordinating by Et is the same as subordinating X by Et and then killing, 
due to the fact that the sample paths of Et are continuous and nondecreas- 
ing. 

COROLLARY 3.4. Let f G C^ k (D) be a 2k-times continuously differen- 
tiable function of compact support in D. If k > 1 + 3d/4, then the equation 
(3.1) has a classical (strong) solution. In particular, if f G C^°(D), then the 
solution of equation (3.1) is in C°°{D). 

Proof. By Example 2.1.8 of [12], \(j) n {x)\ < {\ n ) d/A - Also, from Corol- 
lary 6.2.2 of [13], we have A n ~ n 2 / d . 

Applying the Green's second identity A;-times, we get 



Using the Cauchy-Schwarz inequality and the fact / G C c (D), we get 



{r D (X(E))>t} = {X(E([0,t}))cD} 

= {X([0, E t ]) CD} = {t d {X) > E t } 



(3.18) 






nl/2 



where c& is a constant independent of n. 
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This and equation (3.18) give \f(n)\ < Cfc(A n ) k . 
Since 



A/(x) = J2 -Kf(n)4>n(x), 



71=1 



to get the absolute and uniform convergence of the series defining Af, we 
consider 



n=l n=l 

oo 



< C k ^( n 2/rf)rf/4+l-fc = Cfc ^ n l/2+2/^2fc/ d 
n=l n=l 

which is finite if 2k /d - 1/2 - 2/d > 1, that is, fc > 1 + 3d/ A. □ 

Remark 3.5. Things are rather easier in an interval (0, M), since eigen- 
functions and eigenvalues are explicitly known. Eigenvalues of the Laplacian 
on (0, M) are (nir/M) 2 , and the corresponding eigenfunctions are sm(mrx/M), 
for n = 1, 2, The form of the solution in equation (3.5) on a bounded in- 
terval (0, M) in R was obtained by [1]. 

Let (3 G (0, 1), = (0, oo)xD and define 

Hl(Z?oo) = ju : Dqo — > R : —u(t,x), g^u(t,x), Lu{t,x) G C(I>oo)J, 
H b L {Doo)=n L {D c 



n ^ u: 



oo y 



<g(x)t /3 " 1 ,for some #GL°° (£>),* >oj. 



Next we state and prove an extension of Theorem 3.1 to general uniformly 
elliptic second-order operators. 

Theorem 3.6. Let < a < 1. Let D be a bounded domain with dD G 
C 1,Q and aij G C a (D). Let {X t } be a continuous Markov process corre- 
sponding to the generator L defined in (2.1). Let Tjj(t) be the killed semi- 
group corresponding to the process {Xt} in D. Let {Et} be the process in- 
verse to a stable subordinator of index (3 G (0,1) independent of {Xt}. Let 
f G D{Ld) nC^D) (lC 2 (D) for which the eigenfunction expansion (of Lf) 
with respect to the complete orthonormal basis {ip n '- n > 1} converges uni- 
formly and absolutely. Then the (classical) solution of 

nGHi(D 00 )nC b (D 00 )nC 1 ( J D), 
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(3.19) 



is given by 



(3.20) 



-^pu(t,x) = Lu(t,x), x G D, t > 0, 

u(t,x) = 0, x G dL>, i > 0; 
u(0,x) = f(x), x&D, 



u(t,x) = E x [f(X(E t ))I(r D (X) > E t )\ 
= E x [f(X(E t ))I(r D (X(E))>t)} 

= -/ T D (i)f(x)g /3 (tr 1 ^)r 1 ^- 1 di 



PJo 



T D {(t/lf)f{x)g p (l)dl. 

Proof. Suppose u(t,x) = G(t)F(x) is a solution of (3.19). Substituting 
in the PDE (3.19) leads to 

F(x)^G(t) = G(t)LF(x) 
and now dividing both sides by G(t)F(x), we obtain 

^ G ( t )/ G(t ) = ^M = -,. 



That is, 



(3.21) QjpG(t) = -/iG(t), t>0, 

(3.22) LF(x) = -fiF(x), x£D,F\ dD =0. 

Problem (3.22) is solved by an infinite sequence of pairs (fi n ,^ n ), n > 1, 
where < fi\ < fi2 < /U3 < • • • is a sequence of numbers such that fi n — > 00, as 
n — > 00 , and ip n is a sequence of functions that form a complete orthonormal 
set in L 2 (D) [cf. (2.4)]. In particular, the initial function / regarded as an 
element of L 2 (D) can be represented as 

00 

(3.23) f(x) = £ /(n)Vn(*) 

n=l 

and by the Parseval identity 



(3.24) 



00 

Wfho = E(/») 2 - 

n=l 



FRACTIONAL CAUCHY PROBLEMS ON BOUNDED DOMAINS 21 

Then with \x n determined by (3.22) and recalling that Ep(—^L n t^) solves 
equation (3.21) with fj, = fi n , we obtain 

G(t) = G (n)E (3 (-fi n t f3 ), 

where Go(n) is selected to satisfy the initial condition /. We will show that 

oo 

(3.25) u(t, a:) = £ f(n)Ep{-pL n tP)xl> n {x) 

n=l 

solves the PDE (3.19). 

Define approximate solutions of the form 

N 

(3.26) u N (t,x) = Y / G (n)E f3 (-^)ij n (x), G (n) = f(n). 

n=l 

Step 1. The sequence {it7v(i, OItvgn is a Cauchy sequence in L 2 (D) n 
L°°(D), uniformly in t £ [0, oo). This is proved as in the proof of Theorem 
3.1. Hence, the solution to (3.19) is given formally by 

oo 

(3.27) u(t,x) = Y, f(n)E p (-fi n t^ n (x). 

n=l 

Step 2. As in the proof of Theorem 3.1, it can be shown that the series 
defining u and Lu converge absolutely and uniformly so that we can apply 
the fractional time derivative and uniformly elliptic operator L term by term 
to obtain that u defined by (3.27) is a classical solution to (3.19). 

Step 3. The stochastic representation of the solution as given in (3.20) 
is also valid. We omit the details, as the proof is again similar to that of 
Theorem 3.1. 

Step 4. We have also a decay estimate for u as in Theorem 3.1, namely, 

\H;t)\\2,D<E fi (-XitP)\\f\\ 2 , D . 
Uniqueness is proved as in Theorem 3.1. □ 

Remark 3.7. The uniqueness of the solution to the boundary value 
problem in Theorems 3.1 and 3.6 is shown in [19], Section 4. However, we 
employ a more direct proof of the uniqueness of the solutions. 

Remark 3.8. It may also be possible to prove the existence of L 2 - 
solutions to the fractional Cauchy problem (3.19), as well as the solution 
formula in the last line of (3.20), using Theorem 3.1 in [5]. To this end, we 
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note that Tjj, the killed semigroup on the Banach space L 2 (D), is uniformly 
bounded and strongly continuous, since Tjj is ultracontractive; see [8, 12]. 
However, that result does not yield the form of the stochastic solution (3.20) 
in terms of a killed Markov process, it does not establish classical (strong) 
solutions, and the proof uses very different methods. 

4. Iterated Brownian motion in bounded domains. In the previous sec- 
tion, it is shown that a killed Markov process, modified by an inverse /3-stable 
subordinator, is a stochastic solution to a fractional Cauchy problem on a 
bounded domain in M. d . The inverse stable subordinators with (5 = 1/2 are 
also related to Brownian subordinators [7]. Since Brownian subordinators 
are related to higher-order Cauchy problems, this relationship can also be 
used to connect those higher-order Cauchy problems to their time- fractional 
analogues. In this section, we establish those connections for Cauchy prob- 
lems on bounded domains in M. d . We also establish an equivalence between 
a killed Markov process subordinated to an inverse stable subordinator with 
(3 = 1/2, and the same process subject to a Brownian subordinator. Finally, 
we identify the boundary conditions that make the two formulations identi- 
cal. This solves an open problem in [7]. 

Let Zt = X(\Yt\) be the iterated Brownian motion, = (0, co) x D and 
define 



Recall that ueC k (D) means that for each fixed t > 0, u(t, •) G C k (D). 

Theorem 4.1. Let D be a domain with dD £ C l,a ,0 < a < 1. Let {X t } 
be Brownian motion in and {Yt} be an independent Brownian motion in 
M. Let {Et} be the process inverse to a stable subordinator of index (3=1/2 
independent of {X t }. Let f £ D(A D ) n C^D) n C 2 (D)[c L 2 (D)\ be such 
that the eigenfunction expansion of Af with respect to {4> n :n > 1} converges 
absolutely and uniformly. Then the (classical) solution of 




u e n A 2(D oc ) n c 6 (AxO n c\D), 

d Af(x) 

— u(t, x) = + A 2 u{t, x), x e D, t > 0, 




(4.1) 



u(t, x) = Au(t, x) = 0, t > 0, x e dD, 



u(0,x) = f(x), xeD 
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is given by 

u(t,x)=E x [f(Z t )I(T D (X) > \Y t \)] = E x [f{X(E t ))I(T D (X)>E t )) 

(4.2) 

= E x [f(X(E t ))I(T D (X(E))>t)}=2 / T D (l)f(x)p(t,l)dl, 

Jo 

where Td(Z) is the heat semigroup in D, andp(t,l) is the transition density 
of one- dimensional Brownian motion {Y-t}. 

Proof. Suppose u is a solution to equation (4.1). Taking the n -transform 
of (4.1) and using Green's second identity, we obtain 

(4.3) ^u(t,n) = ^M + X 2 n u(t,n). 



Note that the time derivative commutes with the ^-transform, as 

<g(x)t- l '\ g e L°°(D), t>0. 



g- t u(t,x) 



Taking Laplace transforms on both sides and using the well-known Laplace 
transform formula 

roo +— P 

1 e st dt = s P-l 



/o r(i-/3) 

for (3 < 1, gives us 

(4.4) su(s,n) -u(0,n) = -X n s~ 1/2 f(n) + \ 2 n u(s,n). 

Since u is uniformly continuous on C([0, e] x D), it is also uniformly bounded 
on [0, e] x D. So, we have lim t ^o Jd n (^' x )4>n(x) dx = f(n). Hence, u(0, n) = 
f(n). By collecting the like terms, we obtain 

, A ~ ~ ( , Rn)[l-\ n s- 1 / 2 ] 
(4.5 u{s,n) = -2 • 

For fixed n and for large s, 

f(n)[l - XnS- 1 ' 2 ] _ ^ 1 /2 / - (n)[s l/2 _ An ] 8 -l/2/(„) 



(4.6) u(s,n) 



XI ( s i/2_ An )( s V2 + A n ) ( S V2 + A n )" 



For any fixed n > 1, the two formulae (4.5) and (4.6) are well defined and 
equal for all sufficiently large s. Since 

(4.7) u(t,n) = f(n)E 1/2 (-X n t 1 / 2 ), 

we can see easily from the proof of Theorem 3.1 that u(t,n) is continuous in 
t > for any n > 1. Hence, the uniqueness theorem for Laplace transforms 
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[4], Theorem 1.7.3, shows that, for each n > 1, u(t,n) is the unique contin- 
uous function whose Laplace transform is given by (4.6). Since m u(t,x) 
is an element of L 2 {D) for every t > 0, and two elements of L 2 (D) with 
the same (/> n -transform are equal dx-almost everywhere, we have (4.2) is the 
unique element of L 2 (D) and (4.7) is its ^-transform. 

Observe that for (3 = 1/2, stable subordinator of index 1/2 has the density 
[3], Example 1.3.19, 



01/2 0"0 



V^TTX 3 



Therefore the density pit, I) of \Y t \ is 



(4.8) 2tg l/2 it/l 2 )l 



2/ 



'Airt 



I 2 



As in the proof of Theorem 3.1, we obtain, using (3.16) for (3 = 1/2, taking 
the inverse Laplace transform and then taking the inverse ^-transform, 



(4.9) 



u(t, x) 



n=l 



oo 



(4.10) 



(4.11) 



E/»^l/2(-^ 1/2 )0n^) 

2tg 1/2 (t/l 2 )r 3 dl 



Ln=l 

oo 

n=l 



A "Vn(x) 
A "Vn(x) 



p{t, I) dl 



T D (l)f(x)p(t,l)dl 
E x [f(Z t )I(T D (X) > \Y t \ 



Note equation (4.10) follows from equation (4.8), and then equation (4.11) 
follows by a conditioning argument. 

We know from the proof of Theorem 3.1 and Corollary 3.2 that u given 
in (4.9) can be represented as 



(4.12) 



u(t,x)=E x [f(X(E t ))I(T D (X)>E t )} 

= E x [f(X(E t ))I(T D (X(E))>t)}, 



where Et is the inverse of stable subordinator of index 1/2. 
Now, equations (4.12) and (4.11) lead to 

u(t,x) = E x [f(X(E t ))I{T D {X) > E t )} 

= E x [f(X(E t ))I(T D (X(E)) > t)] = E x [f(Z t )I(r D (X) > \Y t \ 
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As in the proof of Theorem 3.1, we can easily show that the solution u(t, x) 
satisfies all the properties in (4.1), except A 2 u G C(-Doo) and Au G C 1 (D). 
Using (3.13), one can easily show that Au G C l (D). Hence, it suffices to 
show that A 2 u G C{D 00 ). To do this, we need only to show the absolute and 
uniform convergence of the series defining A 2 u. To apply A 2 term by term 
to (4.9), we have to show that the series 

oo 

^/»^(x)A 2 J E 1/2 (-A n t 1 /2 ) 

n=l 

is absolutely and uniformly convergent for t > to > 0. 
Note that the <f> n -transform of A/ is given by 

(f> n (x)Af(x)dx = / f(x)A(f) n (x) 
D JD 



= / f(x)4> n (x) dx = -X n f(n) 

JD 

and using (3.7), we get 

oo oo 

£ i/>)||^)|a„£ 1/2 (- a "* 1/2 ) ^ E l/wll^M c /2 

71=1 (1 = 1 11 

OO 

^ Ct ^ 2 E l/( n )ll^n( X )l A n < OO, 
n=l 

where the last inequality follows from the absolute and uniform convergence 
of the eigenfunction expansion of A/. 

Observe next that the Laplace transform of 



|^(-A^) + ±L - A^-A^) 



for = 1/2 is 



> 1/2 , 1/2 A?a-V2 

1 + A nS ~ 1/2 - -^r— = 0, 



since the Laplace transform of Ei/ 2 (— X n t 1 / 2 ) is g i%^ A t see e 1 ua ti° n (3-3)] 
and the Laplace transform of \j\fxi is s~ 1 / 2 . Hence, by the uniqueness of 
Laplace transforms, we get 

^E 1/2 (-X n t^ 2 ) + - A^v^-Ant 1 / 2 ) = 0. 
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Now applying the time derivative and A 2 to the series in (4.9) term by 
term gives 

d , , Af(x) a2 , , 
— U (t,x)--^U-A 2 u(t,x) 



^i/ 2 (-Ant 1/2 )A 2 ^ n (x) 



2_^f{n) 4> n {x)-^E 1/2 (-\ n t 1 )- 



n=l 
oo 



Tit 



n=l 



— J E 1 / 2 (-A n i 1 / 2 ) + — ^= - A 2 £ , 1 / 2 (-A n t /3 ) 







which shows that the PDE in (4.1) is satisfied. Thus, we conclude that u 
defined by (4.9) is a classical solution to (4.1). 

The uniqueness follows as in the proof of Theorem 3.1. □ 

Remark 4.2. Theorem 4.1 also holds with the version of IBM defined by 
Burdzy [10]. Here, the outer process is the two-sided Brownian motion which 
is defined by X t = Xf for t > and X t = XZ t for t < 0, where Xf and Xf 
are independent Brownian motions started at and Yt is an independent 
one-dimensional Brownian motion. In this case, the IBM takes the form 
Zt = x + X(Yt) and, using a simple conditioning argument, we can show 
that the function 

u(t,x) = E x [f(Z t )I(- TD {X~) <Y t < t d (X+))} 

reduces to equation (4.9) and hence is also a solution to (4.1). 

Remark 4.3. Let / G C 2k {D) be a 2/c-times continuously differentiable 
function of compact support in D. If k > 1 + 3c£/4, then the equation (4.1) 
has a strong solution. In particular, if / G C£°(D), then the classical solution 
of equation (4.1) is in C°°(D). To see this, we use the estimates obtained 
in Corollary 3.4 to get the absolute and uniform convergence of the series 
defining A/: 



n=l 



]T \Mx)\\f(n)\X n <Y.i X n) 1+d/A c{K)- k <cY J {n 2/d ) l+d/A - 

n=l n=l 

and this is finite if k > 1 + 3d/ 4. 
Recall Doo = (0, oo) x D and define 



JU L^eCODoo), Lu G C 1 (l)), 
at 



<g(x)t- 1 / 2 ,g€L oo (D),t>0\. 
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Theorem 4.4. Let D cM. d be a domain with dD G C 1 '", < a < 1, and 
a™ G C a (D). Let {Xt} be the process in M. d corresponding to L defined in 
(2.1), and {Yt} be an independent Brownian motion in R. Let {Et} be the 
process inverse to a stable subordinator of index (3 = 1/2 independent of 
{X t }. Let f G D(L) n C l {D) n C 2 (D)[c L 2 (D)] be such that the eigenfunc- 
tion expansion of Lf with respect to {ip n : n > 1} converges absolutely and 
uniformly. Then, the (classical) solution of 

d Lf(x) 

Tru(t,x) = V- +L 2 u(t,x), x£D,t>0, 
dt K 1 Jrt 



(4.14) 

u{t, x) = Lu(t, x) = 0, t > 0, x G dD, 
u(0, x) = f(x), x G D 

is given by 

u(t,x) = E x [f{X(\Y t \))L{L{T D (X) > \Y t \))]=E x [f{X{E t ))L(T D (X)>E t )} 
= E x [f(X(E t ))L(T D (X(E)) >t)]=2 T D (l)f(x)p(t, I) dl, 



where Toil) is the killed semigroup in D corresponding to L and p(t,l) is 
the transition density of one- dimensional Brownian motion {Yt}. 

Proof. We can show that the solution u is given by the series 

oo 

E/»^/2(-^ 1/2 )^n(x). 
n=l 

Similar to the proofs of Theorems 3.1 and 4.1, we can prove that this series 
can be differentiated term by term to give classical solutions. 

Let U\,U2 be two solutions with initial value / and boundary values 0. 
Then u = u\ — U2 is a solution with initial value equal to 0. 

Use the fact that L is self- adjoint [same as integration by parts since 
Lu,ueC x (D)] to get 

ij) n (x)Lu(t,x) dx = / u(t,x)Lip n (x) dx 
D Jd 

= -fi n u(t,x)ip n (x)dx = -fx n u(t,n). 
Jd 

Using Lu(t,x) = and applying integration by parts twice, we obtain 





^-u(t,n) =nlu(t,n). 
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Since the initial function is zero, we obtain that u(t, n) = for all i > and 
for all n> 1. This implies that u(t,x) = in the sense of L 2 functions due 
to the fact that {ip n :n > 1} forms a complete orthonormal basis of L 2 (D). 
Hence, u(t, x) = for all t > and almost all x € D. Since it is a continuous 
function on D, we have u(t, x) = for all (t, x) G [0, oo) x D, thereby proving 
the uniqueness. □ 

Corollary 4.5. The two types of PDEs in (3.19) and (4-M), with the 
respective initial and boundary conditions, are equivalent. 

Remark 4.6. We discovered in [7] that the two types of PDEs (3.19) 
for (3 = 1/2 and (4.14) are not equivalent without any further assumptions 
on the boundary behavior of the solutions in strict subdomains of We 
get the equivalence of the PDEs (3.19) and (4.14) in bounded domains, if 
we consider the right conditions on boundary behavior of solutions. 

Remark 4.7. It would also be interesting to consider fractional Cauchy 
problems on bounded domains with 1 < (3 < 2, and dtu(t, x) = on x E dD. 
For unbounded domains, some results along these lines can be found in [6, 9]. 
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